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ABSTRACT 

This thesis is a computer simulation of secondary electron 
emission (SEE) from monocrystalline metals and alkali-halides 
using a modified version of the Harrison, Carlston and Magnuson 
single collision theory of SEE. Three cases of SEE are investi- 
gated: the angular dependence of SEE from Cu bombarded by Ar + , 

the dependence of SEE as a function of energy for rare gas ions 
normally incident on the (100) , (110) and (111) faces of 
metal single crystals, and the dependence of SEE as a function 
of energy for Ar + and Ne + ions normally incident on the (100) , 
(110) and (111) faces of KC1. The theory does not accurately 
describe the angular dependence of SEE for monocrystalline Cu 
targets, but does accurately predict the modified sec 0 depen- 
dence found experimentally in polycrystalline studies. For 
the metal targets, the difference between the theoretical 
kinetic secondary emission result and the experimental datum 
is identified as potential secondary emission. The alkali- 
halide SEE simulation agrees reasonably well with experiment. 
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I. INTRODUCTION 



Secondary electron emission is the release of electrons 
from target materials bombarded by primary ions. The first 
observation of secondary electron emission in 1889 was by 
Villard who noticed electrons being ejected from the cathode 
in a discharge tube. Secondary electron mission (SEE) is 
characterized by the coefficient of ion-electron emission, 

Y, which is defined as the number of electrons liberated per 
incident ion. 

Secondary electron emission can be produced by two 
mechanisms, potential emission (PSE) and kinetic emission 
(KSE) . Potential emission arises from the fact that as an 
ion approaches a target material it may become energetically 
feasible for a free electron to be emitted from the target 
surface either by Auger neutralization or Auger recombination. 
On the other hand, kinetic emission arises from the actual 
collision between the electrons of the impinging ion and 
the electrons in the target atom. In general, y depends 
primarily upon the PSE mechanism at ion energies below one 
Kev and the KSE begins to contribute at ion energies above 
one Kev. Although this is a fairly good rule of thumb, 
variables such as the work function of the target material, 
ion charge, surface state of the target, target temperature, 
and whether the monocrystalline or polycrystalline nature 
of the target affect the value of the KSE threshold. 
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Medved et al . [1] have shown that potential emission 



from Ar + and Ar° on Mo increases with energy, with the 
total SEE at higher energies consisting partly of a KSE 
component and partly of a PSE component. In general there 
exists a transition region that is not sharply defined and 
for this reason it is felt that PSE contributes to secondary 
emission in the energy range above one Kev, even though the 
KSE should be the dominant emission mechanism in the one 
to ten Kev energy range. Both types of emission are investi- 
gated in this thesis. 
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II. STUDY OBJECTIVES 



In 1964 Harrison, Carlston and Magnuson [2] proposed 
a single collision theory of secondary electron emission 
from monocrystalline target materials in the one to ten Kev 
primary ion energy range. The original paper develops a 
theory applicable to monocrystalline targets bombarded by 
inert gas ion beams directed normal to the (100) , (110) 
and (111) crystal faces. This thesis will extend this 
theory to explain the following phenomena in the one to ten 
Kev primary ion energy range: 

1. Secondary electron emission from a single crystal 
plane rotated arbitrarily about a crystallographic 
axis, 

2. Secondary electron emission from metal single 
crystals, and 

3. Secondary electron emission from alkali -halide 
monocrystals . 
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III. SECONDARY ELECTRON EMISSION THEORY 



Secondary electron ejection mechanisms are either of 
the potential or kinetic type. Potential ejection can be 
caused by resonance neutralization, resonance ionization, 
Auger neutralization, or Auger de-excitation. Kinetic 
ejection requires the direct transfer of the ion’s kinetic 
energy to the electrons of the target atom. Both mechanisms 
are discussed below. 

Resonance neutralization occurs when it is energetically 
feasible for an electron in the target material to ’’tunnel" 
through its potential barrier and approach the position of 
the incoming ion. This stage of the process increases the 
energy of the ion and target atom system by an amount equal 
to the energy needed to remove an electron from the metal. 
This energy can range from near zero to <f> , the work function 
of the target material. In the second stage of resonance 
neutralization, the electron and the ion combine, effectively 
at infinity, to reduce the system energy to an amouht , 
the ionization potential of the ion. This mechanism is not 
quite as simple as it seems, since it is affected by the 
image potential and coulombic forces of the electronic 
structure of the atoms in the crystal. According to the 
Franck-Condon principle, the electronic transition of the 
electron from the target material to the ion can only occur 
with appreciable probability when the position and velocity 
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of the ion are not affected by the transition. By this 
principle we expect a crossover of the initial and final 
potential energy curves of the system. Taking the image 
potential into account, resonance neutralization can occur 
for reasonable separations of ion and target atom when 



(Vi - 



V - + 



+ I < w 



and 



(V. - V ) + I > <p . 
v 1 e J s Y 

Therefore <p < (V. - V ) + I < W , 

Y K 1 e J s a * 



where <f> is the work function of the target material, , the 
ionization potential of the ion, V , the excitation potential 
of the primary ion, I s , the surface image potential and , 
the initial kinetic energy of the primary ion. 

The second possible potential ejection mechanism is 
resonance ionization which is essentially the reverse process 
of resonance neutralization. In resonance ionization an 
excited atom at infinity is ionized into an electron and an 
ion. The electron then proceeds to the metal surface and 
occupies any of the unfilled available states above the Fermi 
level of the metal. As in resonance neutralization, there is 
a range of ion-target atom separation for which the process 
is most likely to occur. There is also the possibility that 
a highly excited atom will be ionized very close to the 
target surface with a subsequent high probability that a 
secondary electron will be produced. 
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Although the above processes are possible, they are 
not the ones that make the largest contributions to secondary 
electron emission at low energies. Auger processes, which 
occur at distances of approach of 3 X or less, are the 
primary mechanisms contributing to PSE in the low Kev 
energy range. These can be summarized as follows [3] : 



and 



Auger neutralization: 

X + + ne -+ X + e + (n-2)e~ , 
m v J m 9 



Auger de- excitation : 

X m + ne -* X + e + (n-l)e , 
m v J m 9 



where X + and X m refer to the bombarding ion or metastable 
atom, respectively, and ne^ is the number of free electrons 
initially in the target material. In both cases, the final 
state of the system is such that a free electron and a 
neutral atom are produced. Auger neutralization occurs 
when the energy of the initial system, target and ion, is 
greater than the resulting system, X + e + (n-2)e m , at 
infinite separation. Schematically the process is shown in 
Figure 1 [3]. Auger de- excitation , on the other hand, 
involves a metastable atom which approaches the target 
material. As a result, the initial energy curve, shown in 
Figure 2 [3], is not modified by an image potential and 
hence the initial and final energy states never cross. 



v v 
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SCA) 



Figure 1. Potential energy diagram for a helium ion as it 
undergoes Auger neutralization at a tungsten 
surface (3) . 
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5(A) 



Figure 2. Potential energy diagram for a metastable atom (He) 
as it undergoes Auger de-excitation at a tungsten 
surface (3) . 
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These two processes, to a great extent, account for low 
energy (ev range) secondary electron emission. 

Kinetic emission is the result of the direct transfer 
of the incoming ion’s kinetic energy to the atom's electrons. 
Part of this energy is absorbed in the recoil of the atom, 
while some of the remaining energy is transferred to the 
electrons of both the ion and target atom, causing an 
excitation of both of their electronic structures. As the 
two nuclei separate, some of the perturbed electrons are 
liberated as secondary electrons. The ion and recoiling 
target atom are slowed down within the crystal structure 
of the target material and are likely to undergo Auger 
neutralization or de-excitation. If this process occurs 
sufficiently close to the surface of the target material, 
there is a probability that these electrons also will 
escape from the target materials surface and contribute to 
kinetic secondary electron emission. 
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IV. THEORIES CONCERNING SECONDARY ELECTRON EMISSION 



A. THEORY OF PARILIS AND KISHINEVSKII 

Of all the theories proposed to explain kinetic secondary 
electron emission, the theory of Parilis and Kishinevskii 
(PK) seemed to be one of the most promising [4]. PK 
divided the secondary electron emission problem into two 
separate subproblems which were analyzed separately and 
subsequently collated to develop the theory of KSE. These 
two subproblems included: 

1. The search for a sufficiently effective mechanism 
of transfer of the kinetic energy of a moving 
particle to the (target) electron; and 

2. The description of the yield of excited electrons 
from the target. [4] 

PK assumed that the solution to the first problem was of the 
form proposed by 0. B. Firsov for inelastic scattering of 
ion-atom collisons in the Kev range [5]. Secondly, PK 
visualized the secondary electron emission process as one 
in which the primary ion ionizes many atoms as it progresses 
along its trajectory in the crystal. Electrons produced 
in this manner, or by further ionization of the primary 
ion, then proceed to the surface suffering collisions and 
energy losses as well as possible capture in route. The 
combination of the electron production phase and the elec- 
tron escape phase then yields the KSE from the crystal. 

The actual collision between the ion and the atom was 
seen by PK as being accompanied by an overlap of the 
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electron shell structure of the two collision partners, 
during which time there occurred an electronic interaction. 
Energy was transferred between colliding particles, exciting 
their electrons to a higher energy state. To calculate this 
excitation energy, a statistical picture of the electron 
distribution of the particles was used. The theoretical 
relationship, based on the Thomas- Fermi (TF) model of the 
atom, leads to the following formula for the energy trans- 
ferred between the colliding particles: 

2 2 7 

E = — / (/ <J) dS) R • dR , 

4ir z H- i s “ 



where R is the radius vector connecting the centers of the 
colliding particles, S is the surface separating the domain 
of action for the potential of the two particles, and <|> 
is the potential on this surface given by: 



- [- 



CZ-, 



Z 2> e 



] X [ 



1.13CZJ 






TF 



2 2 

where is the charge on each particle, a^p = h /me and 
X is the TF screening function [8]. Physically, the inter- 
action can be thought of as analogous to the collision of 
two tennis balls. As they impact, both deform, with the 
vector joining their centers being perpendicular to the 
area S of the plane of interaction which separates them. 
Obviously the total energy transferred is a function of the 
relative velocity with which the ion and atom impact and 
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the electron distribution associated with each. In terms 



of the tennis ball analogy, the greater the relative impact 
velocity, the larger the deformed area S, the shorter the 
radius vector distance R, and the greater the number of 
electrons involved. 

The analysis was rewritten by PK in terms of the impact 
parameter associated with the collision. Hence by changing 
variables 



where u is the initial velocity of the ion, E is the 
excitation energy of the colliding particles, p is the impact 
parameter and V(r) is the potential of repulsion between the 
particles . 

To obtain the correct behavior of the emission curve, 
the cross section of electron ejection into the conduction v 

band had to be calculated. Because more than one electron 
may be lifted from the valence band to the conduction band, 
the cross section was determined to be 







E(P) 




X 2 (p) dp 



P 



o - Z (cons tant) • (Probability of ionization) 



P 



•(Impact parameter) 
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which in the integral limit approaches 



a 



2tt 



r P l E(p) 

i J 



P dp . 



Here is the impact parameter at which E(P^) = 6 - , 

where 6 is the energy depth of the filled conduction band, 
( p is the work function of the target material and J is the 
average ionization potential for the outer shells of the 
atom. For 1/4 < < 4 PK found the cross section to 

be given by 



a(u Q ) = 



1 . 39 a T p h 



( Z 1 ♦ z 2 ) 



(Z 1 1/2 + z. 



\rm s ^ u o ) 



and S(u q ) is approximated by 



S(u o ) 



5.25 u arctan 
o 



(.6 x 10' 7 (u 



- u . ) ) 
mm' ' 



where u . is the threshold velocity. After determining 
the probability that an electron will be promoted from a 
filled band to the conduction band, PK assumed that most 
of these electrons remained in the conduction band because 
of the surface potential of the metal target. However, 
because a hole is left in the valence band, an Auger recom- 
bination is possible with subsequent emission of an electron 
from the metal. The probability that this process occurs is 
given by PK as 

w (6 ) = 0.016 (6 - 2<J>) . 
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PK assumed that both PSE and KSE were due to Auger 
processes [4] . 

After the secondary electrons have been produced, they 
suffer collisions as they proceed to the surface of the 
metal. These collisions degrade their energy and increase 
the probability that they will be lost through recombination. 
PK assumed that the number of electrons which reach the 
surface is given by an exponential law of the form exp (-x/A) 
where x is the depth at which the electrons are formed in 
the crystal and A is the mean free path of the electrons 
in the material. Specifically, the number of electrons 
emitted per ion incident on the surface is given by 

Y = £ (Number of atoms , (Prob. that an electron 
x per cubic centi- will be ejected into 

meter) the conduction band 

and a hole left in 
valence band) 

•(Prob. that an electron 
will escape to the 
surface of the crystal) 

which in the limit goes to 

X 

Y = f n N a (u) w (6) exp (-x/A) dx , 
o 

where X is the maximum depth at which the ion retains 
n 

sufficient energy to ionize a target atom. Upon taking 
into account the degradation of energy of the primary ion 
as it suffers multiple collisions and the diffusion cross 
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section of the electrons, PK arrived at the final expression 
for y : 



Y = N w (6) X [o (u ) - Aa (u ) ] , 
o o 

where 

Aa (u Q ) = exp (-u /(KA) / U o exp (u 2 /(KA)) du 

u min 

and 

2.487T N • a^ c e 2 Z-. Z 0 
^ _ TF 1 2 

(M x ♦ M 2 ) (Z ^ 17 2 + Z 2 1/2 ) 2/3 

The PK theory predicts the following trends [3] : 

"1. Low velocity region: u = u . , Aa (u ) term is 

J 6 o min’ K o J 

important..., y increases slowly with u , 

(a (u ) - Aa(u )) = u 2 - (3/2 u . ) 2 . Hence 
vv o v o o v mm 

? 

initially Y a u o , i.e., y is a linear function 
of energy . . . 

2. High velocity region: (g(u q ) - Aa(u Q )) = 

- 1 - 7 

u tan (.6x10 (u -u . )). This expression 
o v K o mm" r 

asymptotically approaches the straight line 

v 

a (u ) - c (u -u . ) and therefore increases 
v o o mm 

linearly with velocity. Extrapolating the 
linear portion back to a(u Q )=0 a value of the 

7 

intercept u^=l. 05x10 cm/sec is obtained; 
this is independent of ion target combination. 

3. Very high velocities: Here ion penetration depths 

are greater so that electrons are formed further 
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from the surface. The number escaping to the 
surface and therefore the yield is expected 
to pass through a maximum. This is confirmed 
experimentally . 

4. No yield dependence on ion charge is predicted. 

5. A dependence of yield on the ion target combination 

of the form {{ 1 ^+ 1 ^) / { 1 ^ ^ + 1 ^ for heavy 

ions and ( + for light ions. 

If the term Aa(u Q ) is important ... the yield is 
independent of ion-target combination to a first 
approximation. 

6. Dependence of yield on the angle of incidence 

G is expected to be of the form y^ = sec 0 since 
the probability of electron escape is a function of 
the shortest distance to the surface and the proba- 
bility of formation of electrons is a function of 
the actual distance traversed. 

7. The effect of the isotopic mass on yield is 
anticipated to be connected with the different 
retardation rates for the isotopic pairs. n [3] 

All of these conclusions are approximately confirmed by 
the experimental results. 

B. THE THEORY OF HARRISON, CARLSTON AND MAGNUSON 

The PK theory uses many simplifying assumptions which 
introduce approximations into the results obtained. In 
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particular, the PK theory predicts secondary electrons of 
only average energy, two to five ev, and data obtained 
by Wolff [6] of SEE from copper bombarded by ten Kev argon 
ions indicates some secondary electrons emitted have energies 
greater than 140' ev. 

In the PK theory secondary electron emission was viewed 
as a two-step process in which ions were first produced 
along a path in the target material and subsequently diffused 
to the surface with possible scattering or Auger emission 
occurring along the path. This model does not yield the 
high-energy secondary electrons which are reported experi- 
mentally. In an effort to explain this discrepancy between 
experimental and theoretical result, Harrison, Carlston 
and Magnuson [2] (HCM) proposed a ’’single collision" theory 
of secondary emission. This theory, though in many respects 
similar to the PK theory, has as its basic philosophy an 
electron production mechanism which differs from the model 
proposed above. 

HCM proposed a model in which an ion approaching a target 
surface first undergoes resonance neutralization and subse- 
quently collides with the first layer of atoms in the target 
material in the un-ionized state. Even though the primary 
particle is neutralized before the collision, the terminology 
of "ion" is retained here to distinguish between the primary 
particle and the target atom after collision. For simplicity 
HCM assumed that the electrons liberated due to the inelastic 
energy transferred between the electronic structures of the 
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ion and the atom all come from the ion, which allowed them 
to ignore the energy partition between collision partners. 

Theoretical work by narrower [7] indicates that the Auger 

- 14 

process requires a time on the order of 10 seconds while 

in the one to ten Kev energy range the collision process 

- 1 4 

takes on the order of 5 x 10 seconds. This indicates 
that, statistically, the ion usually does not have time to 
neutralize before its next collision and for all practical 
purposes cannot contribute further to the secondary electron 
emission process. These results lend credence to the single 
collision model. Mathematically HCM proposed: 



(hk£) 

Y 

KSE 



(hkfc) 

where y is 

KSE 



(hk£) 

V , S MAX „ , p , D 

K / n e (s , Ejpp) P 



(hU) 



(s) ds , 



the number of electrons emitted per incident 



ion when the ions are normally incident upon the (hk£) surface 
s is the impact parameter of the ion on the atom; E,ppp is 
the inelastic energy transfer between the ion and the atom; 
n r s n ) i s the number of electrons this energy will 

produce; and is the probability that the impact 

parameter s will occur in the (hkJl) surface. This equation 
contains three adjustable parameters, the constant K and two 
constants which enter into the calculation of 4, the electron 
density, and the interaction potential V(r). By symmetry, 
these last two constants must be the same for all orientations 
of the particular crystal plane, (hk£) , in question. 
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According to HCM, the constant K is a factor which accounts 
for three physical processes: 

1. The probability that an electron will be emitted 
from the surface (a value close to 1/2); 

2. The probability of additional emission by an 
Auger process; and 

3. The probability that a high energy electron 
escapes before it is scattered back into a 
conduction band. [2] 

All of these processes are orientation independent. The 
factors which appear in the formulation of the HCM secondary 
electron emission equation are discussed in detail below. 

1 . The Inelastic Energy Transfer 

In the HCM theory Eppp is the inelastic energy 
transferred in the collision process and is basically the 
same as E(p) in the PK theory. However, as HCM point out, 
because the interaction involves atoms spaced regularly in 
a lattice array, the inelastic energy integral does not have 
an infinite upper limit as was assumed by PK. Specifically, 



e tff( s > e 0 ) “ K tfF 



max 

f (l-V(r)/E o ) dr 

J (l-V(r)/H o -s 2 /r i ) 



Tn 



# 



r EM 
r/2 



1, 2 (P) dp 
P 
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where 



V - 

TFF 



-K (E /(2m )) 1/2 (Z.+Z ) 2 
lra H 



a H = h 2 / (me 2 ) = . 529 X 



and E q is the energy of the ion in the center of mass system 
Z i and are the ion and atom nuclear charges, respectively 
r min t ^ ie distance °f closest approach; and <j> is the 
Thomas-Fermi (TF) screening function: 



w ^ _ ri + fo n >>0.8034.-3. 734 
4 > Cp J - [l + ( a pp) ] 



where 



a F = 



a Tp ((Z 1 +Z 2 )/144) 



1/3 



0.8853 



and a,pp is one of the adjustable constants introduced above. 
By experimenting with different sets of potentials, in 
conjunction with the screened electron distribution, HCM 
found that the best agreement with experiment was obtained 
when the potential V(r) was of the Born-Mayer form [2]: 

V(r) = exp (A+Br) = A* exp (Br) 

with B always negative, and where A is the second of the 
above mentioned arbitrary constants. HCM referred to the 
integral containing dp as the electron integral. The con- 
vention is retained here. 
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Several important results were noted by HCM in the 

formulation of their original model. First, the dynamical 

f si 

integral was extremely sensitive to its upper limit r^ a ^. , 
whereas the electron integral was not sensitive to varia- 
tions in r C w. For this reason r„., was set at a value some- 
what larger than the nearest neighbor distance R q , in the 
crystal lattice, and retained this value for all orienta- 
tions of the crystal lattice. On the other hand, the upper 
limit of the dynamical integral was extremely sensitive to 
the nearest neighbor distance. Such a behavior is expected 
because as each ion collision with a particular atom termi- 
nates, another must begin with one of its neighbors. From 
a particle dynamics point of view, HCM found that by 
treating the actual collision with the lattice as a binary 
collision between the incoming ion and any single atom 
within the first repeat distance of the crystal plane, the 
theory yielded results which approximated those obtained by 
experiment. Figure 3 indicates the geometrical relation- 
ship between various quantities associated with the 
dynamical integral [2]. 

2 . Collision Parameters 

Because of the repetitive nature of the lattice, or 
lattice plane, it is not unreasonable to think of a "repre- 
sentative area" which characterizes the lattice plane as a 
whole and hence its distribution of impact parameters. If 
no lattice were present, and only one target atom, there 
would be no characteristic area and the occurrence of a 
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Figure 3. Relationships between various quantities appearing 
in the Dynamical Integral. 
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collision with any impact parameter would be possible, 
yielding an infinite collision cross section. As an ion 
approaches a plane of atoms in a crystal, however, 
there exists .an allowable set of impact parameters for the 
ion-atom collision. Specifically, in a lattice there exists 

fhk£! 

an impact probability distribution, p v J (s) , for each 
crystal orientation. For normal incidence HCM proposed the 
areas shown in Figure 4. These areas are representative of 
the normal impact of an ion with any atom on the specified 
crystal surface, viz., by symmetry operations they fill plane 
space and are representative under any of these symmetry 
operations of what an ion "sees" as it approaches the planar 
surface normally. There are many distributions which could 
be considered probability distributions. HCM considered the 
following : 

1. The distribution of values of s measured to the 
nearest surface atom; 

2. The distribution of all values measured to all 
atoms which project their area into the repre- 
sentative area; 

3. The distribution of smallest impact parameters 
when all atoms are considered [2], 

The third case was the only one which yielded results in 

agreement with the experimental data. The distributions 

for normal incidence on FCC and BCC metal targets are shown 

in Figures 5 and 6. 

Operationally, the distributions are determined as 
follows: First, the representative area is chosen. 

’’Points" are shot at the representative area, and the distances 



32 




Figure 4. Representative areas for normal incidence. 
After (2) . 
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NORMALIZED PROBABILITY DENSITY 




(Ro units) 



Figure 5. Distribution of impact parameters for normal 
incidence on face -centered cubic crystals. 
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NORMALIZED PROBABILITY DENSITY 




(Rounits) 



Figure 6. Distribution of impact parameters for normal 
• incidence on body-centered cubic crystals. 
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between the impact point and the centers of the atoms lying 
in the representative area are measured. For each atom the 
distances for points that are measured to be closest to it 
are tallied. After a sufficient number of points are tallied 
to yield good statistics, the distribution for each atom is 
weighted by a factor to account for the fractional portion 
of atom lying in the representative area. The atoms' dis- 
tributions are then added and scaled so that the total 
area under the resulting curve is unity. 

3 . Number of Electrons Emitted- - Russek Model 

One of the major shortcomings of the PK theory was 
that it used average ionization energies to obtain an ioni- 
zation cross section. This is to say the PK theory did not 
consider the electronic structure of the target atom. HCM 

adopted a model proposed by Russek [9] to obtain the number 

(hkjn 

of electrons liberated per collision, n^ ^(s,E,ppp). In 
his original paper, Russek proposed two ionization models of 
interest, the uniform ionization model (UI), and the staggered 
ionization model (SI). The UI model assumed that as electrons 
were liberated by the collision process, the amount of energy 
required to liberate any one electron was the same as that 
required to liberate any other electron. On the other hand, 
the SI model assumed that each electron liberated required 
more energy than was required to liberate the preceding 
electron. For the SI case Russek shows that the probability 
that n electrons will be liberated from a shell, which 
contains M electrons is given by: 
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1-nE 



ION 





n 



where k £ (E 




nE^^/E^ON) < ^+2 and M is the number 
n+1 n+1 “ 



of electrons available in the outermost shell of the atom; 



electrons and E T pp is the energy available to cause this 
removal. Two things should be noted. First, the impact 
parameter enters this calculation through the inelastic 



then, from probability arguments, the number of electrons 
1 iberated is : 



HCM were unable to obtain good results with the Russek 
UI model. In all cases, y increased much too rapidly as a 
function of energy. On the other hand, the SI model yielded 
reasonably good results, where spectroscopic values for E*^ 
were used. 



E^^ is the energy required to remove the n^ 1 of these 



n 



energy transfer. Second, when E"^^ = E^^ one obtains 

n+1 n 

the UI model. After normalizing the function, , viz. 





\ 
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4. HCM Results 



HCM investigated the Ar - Cu system in detail. 

They first examined the Ejpp dynamical integral for various 
types of potential functions and found that the PK potential 
which was of the. Thomas - Fermi- Firsov type was too hard for 
small values of E , viz., the hard core radius of the target 
atom was too large, in this energy range, to yield accepta- 
ble results. Attempts at softening the potential by decreas- 
ing a^p failed. For this reason HCM assumed a potential 
function of the Born-Mayer type. The potential function 
used in each case was found by matching the exponential to 
the TF screening function at some separation. 

HCM assumed that the experimental data contained both 
potential, Yp S p> and kinetic, Y K g E > contributions to the 
total secondary electron emission coefficient. Therefore, 

HCM identified the y at which KSE begins to contribute as 
Ypsp and assumed this value remains constant for all energies. 
This value for Ypgp was then subtracted from the experi- 
mentally determined y for the (111) surface at 10 Kev. The 
(111) value of y^gp at 10 Kev, obtained from theory, was 
scaled to this point, and all other points were multiplied 
by the scale factor used in performing this operation. 

Finally the constant value of Ypgp was added to all the 
scaled theoretical values. The result is the HCM value for 
the total secondary electron emission. The results of the 
HCM theory for the Ar + - Cu system are presented in Figure 7. 
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SEE COEFFICIENT 




Ar + - Cu results obtained from the HCM original 
model. Taken from (2). 



Figure 7. 



In general, HCM's theory agrees reasonably well with 
experimental results for all systems studied except Ar - Mo. 
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V. MODIFICATIONS TO THE HCM THEORETICAL MODEL 



A. INTERATOMIC POTENTIAL 

Independent of the philosophy behind the basic interaction, 
the HCM theory, like the PK theory, is highly dependent on two 
specific assumptions. First, both assume that the inelastic 
energy transferred is of the form proposed by Firsov. Second, 
both theories must assume an interatomic potential function 
which is consistent with the electron distributions of the 
ion and atom. Since the Firsov model is presently the only 
comprehensive model available for the inelastic energy 
transfer, it has been retained here as an approximate model 
of the inelastic collision process and is assumed to be 
sufficiently accurate for our purposes. However, the inter- 
action potential used in the Firsov model need not be of 
the Thomas- Fermi form, but may be derived from a more 
accurate representation of the electron distributions of the 
ion and atom. 

The interaction potential is derived from the electron 
distributions of the collision partners. These potentials 
exist on three general levels of sophistication: (1) The 

Bohr and Born-Mayer potentials which are derived from 
electron distributions which provide screening of the 
nucleus in a self-consistent, but not necessarily accurate, 
manner; (2) Thomas-Fermi and Thomas - Fermi - Dirac potentials 
which are derived from statistical pictures of the electron 
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distributions of the collision partners; and 

(3) Potentials generated from electron densities derived 

numerically from Schroedinger ' s equation. 

Parilis and Kishinevskii used a Thomas-Fermi potential 
function in the -theory they proposed. Statistical models 
of the atomic electron distribution, such as in the TF 
model, are based on studies of inhomogeneous gases and 
approximate solutions of Schroedinger ’ s equation consistent 
with this picture [10, 11]. They require that the total 
electron energy be a minimum and that the electron density 
integrate to the total number of electrons in the system. 
Approximations may be made to evaluate the exchange and 
correlation energies. When the exchange and correlation 
energies are omitted and the kinetic energy is approximated 
by a plane wave, the TF atomic model results [11]. If the 
exchange energy is retained and the kinetic energy is 
approximated by a plane wave, the TFD atomic model results 
[11]. The TF model of the atom has been used extensively 
by Firsov, whereas the TFD model appears in papers by 
Abrahmson [13]. 

Potential functions based on Hartree-Fock solutions of 
the atomic Schroedinger equation have been discussed by 
Wedepohl [12] , Harrison [13] , and Wilson [14] . These 
potentials assume the adiabatic approximation, use 
spherically symmetric electron distributions which have 
shell structure, and take into account the correlation, 
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kinetic and exchange energies through various approximations. 
They appear to lead to improved calculations of the interac- 
tion potential [8, 14, 15]. The present investigation 
utilizes an interaction energy calculation similar to that 
proposed by Wedepohl [12]. However, the electron densities 
from which the potential is derived are calculated from the 
Hartree- Fock-Slater self-consistent wave functions calcu- 
lated from the Herman- Skillman [16] computer program. The 
reader is referred to the original book [16] for detailed 
information on the calculation of the atomic electron dis- 
tribution functions and to the work by Torrens [8] for a 
general discussion of the interaction potential. 

B. THE BASIC INTERACTION 

Harrison, Carlston and Magnuson assumed that all of the 
secondary electrons produced in the single collision model 
came from the ion. This assumption was tested by choosing 
various values of E*^ N in the Russek model and noting the 
dependence of n e ( E ppp) on E J° N - This experiment was car- 
ried out for the Ar + - Cu system. The results were most 
favorable for E* 0N equal to 15.8 ev, the first ionization 
potention of argon, and least favorable for neutral copper’s 
first ionization potential. This approximation was an 
attempt to simplify the model so that detailed models of 
energy sharing mechanisms could be avoided. In this 
present investigation, this assumption has been relaxed to 
include contributions to the secondary emission by the 
target atom. 



43 



Theoretically, the interaction is viewed as the collision 
of two atoms having spherically-symmetric , shell structured 
electron distributions. As the atoms collide, they are 
assumed to form a "molecular" structure whose electron 
ionization potential distribution is comprised of the 
separate ionization potentials of the electrons of the atoms 
involved in the collision. Since electrons are ejected on a 
"least energy" principle, the values of E*^ used in the 
Russek model have been modified by include both primary ion 
and target atom ionization energies which are arranged in 
increasing order of magnitude, irrespective of the atom to 
which they belong. The mean number of electrons is a compli- 
cated function of the inelastic energy transfer and the 
energy necessary to effect each electron's removal. 

C. THE DISTRIBUTION OF IMPACT PARAMETERS 

To accomplish the objectives of this present investigation, 
one additional modification to HCM’s original paper had to be 
made. Since the angular dependence of y was to be studied, 
a new distribution of impact parameters had to be calculated. 
Fortunately, according to the HCM paper, this modification 
to the mathematical model is the only change required to 
take into account rotations of the crystal plane around an 
arbitrary crystallographic direction. No other factor 
entering into the secondary electron emission equation is 
orientation dependent [2]. 
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Two possible solutions to this problem are proposed 
below. First, a crystal plane normal to the direction of 
incidence of the ion beam can be found, from which a repre- 
sentative area can be constructed. Such a formulation would 
be analogous to cutting the crystal in such a manner that 
it presented a planar area normal to the incident ion beam 
and finding an area on this plane which would be representa- 
tive of that orientation. Alternatively, a representative 
area constructed from the crystal face which was experi- 
mentally rotated can be used and mathematically rotated 
about the specified crystallographic axis. In practice, the 
first method is computationally restrictive because it 
requires that a representative area on a plane normal to the 
incident ion beam be found for each orientation. For this 
reason it is not pursued further here. 

The second method mentioned above yields a result 
that is more in line with the reasoning behind HCM's 
original definition of the distribution of impact parameters. 
Practically speaking, a representative area is found for 
each crystal face which fills plane space and is representa- 
tive of what an ion sees as it approaches the crystal plane 
from an arbitrary direction. Unlike the representative 
areas for normal incidence, angularly dependent representa- 
tive areas must be truly representative for all possible 
rotations of the plane and not simply for one particular 
orientation. The representative area for the FCC (100) 
surface which can be rotated in either a positive or negative 



45 



sense about the < 1 1 0 > crystallographic axis is shown 
in Figure 8. Note that atoms which lie within the first 
repeat distance of the crystal plane are eligible target 
atoms in the representative area by HCM’s definition of 
representative area. Also note that these atoms from the 
lower atomic plane appear to move as the crystal plane is 
rotated in a positive sense about the <110> axis (see 
Figure 9) . This behavior can be verified by physically 
rotating a stick and ball model of the lattice structure 
about the indicated axis. 

Once the representative area is chosen, its rotational 
properties and corresponding distribution of impact 
parameters are calculated as follows. The representative 
area is placed in the first quadrant of a rectangular 
coordinate system (x,y) (see Figure 10). In this system, 
the crystallographic axis of rotation is specified by two 
parameters : 

1. The angle of the crystallographic axis makes with 
the x-axis of the rectangular coordinate system 
(x,y); and 

2. The y-axis intercept of the crystallographic axis 
in the rectangular coordinate system (x,y). 

This system of coordinates is then transformed into a primed 
rectangular reference frame (x’,y f ) for which the crystal- 
lographic axis of rotation is the x’-axis, and whose origin 
is the point at which the crystallographic axis intercepts 
the y-axis, (* R >y R ) ( see Figure 10). The equations relating 
the transformation of points from the (x,y) coordinate 
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Figure 8. Representative area of a (100) FCC crystal face 
rotated about <110> in either a positive or 
negative sense. 
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Figure 9. Representative area of the (100) face of an FCC 
crystal rotated in a positive sense about the 
5 < 1 1 0 > direction. 
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Figure 10. Rotation geometry for a (100) face of an FCC 
crystal about an arbitrary crystallographic 
direction . 
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system to the (x',y') coordinate system are: 



x’ 


- R cos 


(9 


* <10 - (y R sin 


<f> + x R cos <f>) , 


y 


= R sin 


(0 


- <t>) + (x R s in 


<P - y R cos <|>) . 



where R = (x 2 + y 2 ) 1 / 2 

G = arctan (y/x) 

and <f> is the angle that the crystallographic axis makes 
with the x-axis. If the system is rotated through an angle 
a about the x'-axis, distances to any point (x',y') appear 
to be shortened by an amount cos a in the y' dimension, 
leaving the x' dimension unchanged. Also, if the atom whose 
projected position lies at (x',y') before rotation is an 
atom whose actual position is units below the plane being 
rotated, then its new apparent position on a plane perpen- 
dicular to the direction of incidence after rotation is 
given by the doordinates (x^y 1 cos a - R & sin a). In 
general then the final transformation equations become: 

x’ = R cos' (0 - <f>) - (y R sin <f> + x R cos cf>) , 

y* = (R sin (0 - <f>) + (x n sin <f> - y D cos 0)) cos a - R sin a . 

k k a 

The geometry implicit in these calculations is shown in 
Figure 10. 

The actual calculation of the distribution of impact 
parameters is performed as follows. A mathematically con- 
structed, rectangular grid is placed over the area and 
points which lie within the representative area are 
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systematically chosen. Using the transformation equations 
and the known position of the atoms in the representative 
area, the distance from the points chosen in the area to 
each respective atom are calculated in (x’,y') coordinates 
and the shortest' distance is chosen. The shortest impact 
parameter is tallied against the atom with which it is 
associated, in an array. In the case of equal distances 
between two atoms, the atom with which the ion impacts is 
chosen by a uniformly distributed pseudorandom process [17]. 
The final array contains the distribution of all impact 
parameters tallied with each atom in the representative 
area. To obtain the final distribution of impact parameters 
for the crystal face, the individual distributions are 
weighted, added, and the area under the curve is normalized 
to unity. The number of points used in calculating these 
distributions are on the order of 100,000 per representative 
area. The final unnormalized distribution of impact parame- 
ters for (100) FCC crystals for various angles of rotation 
is shown in Figure 11. Note that all dimensions used in 
calculating these distributions are in terms of R , the 
nearest neighbor distance in the crystal lattice. 
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Figure 11. Unnormalized distribution of impact parameters 
for FCC(IOO) rotated through various angles 
about <110>. 
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VI. NUMERICAL METHODS AND COMPUTER PROGRAMS 



The theory presented in the preceding chapters has been 
incorporated into five FORTRAN IV computer programs which 
are presented in the appendix to this thesis. They appear 
in the order of their use. Specifically, these programs 
are : 

1. The Herman-Skillman radial -electron density program, 

2. The Harrison potential program, 

3. Distribution of impact parameters program, 

4. Numerical interpolation program, and 

5. The Harrison KSE program. 

To calculate the coefficient of ion-electron emission 
(SEE coefficient) , these programs are utilized in the fol- 
lowing manner. The ion- atom- radial -electron densities are 
calculated using the Herman-Skillman computer program [16]. 
These densities are put into the Harrison potential program 
which calculates the interatomic potential, V(r), used in 
the calculation of the inelastic energy transfer, Eppp. 

The distribution of impact parameters program calculates 
the distribution of impact parameters for each crystal face 
under investigation.- Since this is a statistical process, 
the distributions recovered from this program are "smoothed" 
using a first and tenth degree Legendre orthogonal least 
squares interpolating polynomial, which is generated by 
the numerical interpolation program [18]. Finally the 
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smoothed distribution of impact parameters and potential 
function are inserted into the Harrison KSE program which 
has been modified to include all the changes to the original 
theory presented in the "Modifications” chapter of this 
thesis. This program performs the actual calculation of 
the SEE coefficient using the HCM model. The input data 
format to each of these programs is shown in Figures 12 
through 16. The definitions of the input variables used 
in the programs are presented in Appendix A to this thesis. 
The results presented in the following chapter were run on 
the Naval Postgraduate School IBM/360 Computer System. 
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Data Input for the 
Electron Density Program 




Figure 12. Deck setup for running the Electron Density 
Program. 
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Data Input for the 
Potential Program 




Figure 13. Deck setup for running the Potential Program 
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A5(I» J) 



SCALEF,RNA,RCL 



ICODE > A , B , INUMBR 



CHAR 



Data Input for the 
Impact Parameters 
Program 



Figure 14. Deck setup for running the Distribution of 
Impact Parameters Program. 
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Data Input for the 



Interpolation 

Program 



Figure 15. Deck setup for running the Interpolation Program. 



Data Input for the 
KSE Program 



ZB ,ZT , BMAS ,TMAS ,EOS , KEOM, AO, LSKIP , 
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EA(I, I = 1,10) 



EA( I, I = 11,18) 



OUTPUT FROM POTENTIAL PROGRAM 
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HEADER CARD 




FLP(I). . .SMOOTHED DISTRIBUTION OF IMPACT PARAMETERS FOR 
( 100 ), ( 110 ), ( 111 ) 

LCRYST , LTPE, BMAX(I) ,PUD( I ) , ARO( I ) ,LM(I ) , KM,FACL 



Figure 16. Deck setup for running the KSE Program. 
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VII. RESULTS 



Three separate cases of secondary electron emission 
were investigated using the modified HCM model. These were: 

1. Angularly dependent secondary electron emission 
from the (100) face of a copper single crystal 
rotated about a <110> axis; 

2. Secondary electron emission from Ar, Ne , Kr and 
Xe ions normally incident on the (100) , (110) 
and (111) faces of Cu, Ag , Mo, and A1 single 
crystals; and 

3. Secondary electron emission from Ar and Ne ions 

normally incident on the (100) , (110) and (111) 

faces of KC1. 

In the HCM theoretical model, the primary ion is assumed 
to be neutralized before it strikes the target surface. The 
target atoms in the metal are assumed to be ion cores in a 
sea of valence electrons. For these reasons, the actual 
results presented here are obtained by modeling the binary 
collision process as a collision between a neutral, inert 
gas atom and an ion core of the target atom. In the cases 
of copper and silver, the ion cores are fairly well defined, 
since both appear in column I.B. of the periodic table. 
However, in the cases of molybdenum and aluminum which 
appear in the VI. B. and III. A. columns of the periodic 
table, respectively, there are several possible ion core 
configurations that could be used in the HCM model. For 
this reason, various degrees of ionization of the molybdenum 
and aluminum target atoms were used in computing the SEE 
coefficient for these target atoms by way of the HCM model. 
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In the cases of copper and silver, results indicate 
that the best ion core configuration for both target 
materials is the singly ionized state. Note that all graph 
titles refer to the HCM model used, whereas the figure 
captions indicate the physical system actually investigated. 

A. ANGULARLY DEPENDENT RESULTS 

The results obtained from the bombardment of Cu single 
crystals by Ar ions for various angles of rotation about 
the < 1 1 0 > crystallographic axis are shown in Figure 17. 

It is interesting that the HCM model did not predict the 
correct behavior of the SEE coefficient for monocrystalline, 
angularly rotated targets, but did accurately reproduce the 
results obtained for polycrystalline targets [19] . This 
seems to indicate that the single collision model is not an 
accurate one for arbitrary rotations about the crystallo- 
graphic axis, but needs to be modified to account for a more 
complicated interaction with the crystal lattice. 

B. METAL TARGET RESULTS 

An attempt was made in this thesis to verify those 
calculations performed by HCM in their original paper. Not 
too surprisingly , the results obtained here differ appre- 
ciably from those obtained by HCM using their original 
model. Each system is discussed in detail below. 

1 . Copper Bombarded by Ar, Ne , Kr and Xe Ions 

+ 

In HCM’s original paper, the Ar - Cu system was 
treated in detail. Using their original model, they 
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ANGLE OF ROTATION (Deg.) 



Figure 17. SEE coefficient vs. angle of rotation for (100) 
FCC Copper crystal. Direction of rotation was 
* about <110>. 



obtained the results shown in Figure 7, which show excellent 
agreement with experiment. After altering their model, 
using the modifications presented in Chapter IV, the results 
do not agree so well with experiment as did those obtained 
by HCM. (See Figures 18 through 21.) 

In an attempt to explain this discrepancy, the 
difference between the experimentally observed SEE coefficient 
and that obtained here was plotted as a function of energy. 
Theoretically, if the results obtained in this present in- 
vestigation are to be physically meaningful, this difference 
must be due to the potential ejection mechanism operating at 
higher energies. Experimentally, Medvel e t al . [1] have 

ascertained that for Ar bombardment of polycrystalline Mo 
targets, the potential ejection mechanism is still operating 
at energies up to 2.5 Kev, and the PSE contribution is still 
tending to increase with increasing energy. If this is in 
fact the case, the results obtained here have physical 
significance. (See Figures 22 through 25.) 

It is interesting to note that the PSE contribution 
from all three faces investigated exhibit similar behavior 
in that they all peak at approximately the same energy, and 
tend to decrease at higher energies. The exact location of 
the peak, however, varies with the primary ion type. Further- 
more, in all four primary ion cases, the PSE contribution 
curves for the (100) and (110) faces cross at approximately 
seven to nine Kev, indicating that the PSE mechanism may 
be somewhat independent of ion type, while showing some 
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Figure 18. Cu bombarded by Ar + ions normally incident on 
the (100), (110) and (111) planes. 
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Figure 19. Cu bombarded by Ne ions normally incident on 
the (100), (110) and (111) faces. 
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Figure 20. Cu bombarded by Kr ions normally incident on 
the (100), (110) and (111) faces. 
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Figure 21. Cu bombarded by Xe ions normally incident on 
the (100), (110) and (111) faces. 



67 




1 2 3456789 10 
ENERGY (Kev) 



Figure 22. PSE contribution to secondary electron emission 
for the Ar - Cu system. 
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Figure 23. PSE contribution to secondary electron emission 
for the Ne - Cu system. 
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Figure 24. PSE contribution to secondary electron emission 
for the Kr - Cu system. 
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Figure 25. PSE contribution to secondary electron emission 
for the Xe - Cu system. 

«» r. 
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dependence on the plane being bombarded. Such a crossing 
behavior in the PSE contribution as a function of energy 
would not be expected simply by analyzing the SEE coefficient 
curves, since these curves do not cross at any energy. If 
this is the case, then it would follow that the theoretical 
results obtained for each crystal face should be indepen- 
dently fitted to the experimental data, without regard to 
the other planes being bombarded. In most cases this would 
yield results that are in closer agreement with experiment 
and yet do not eliminate the PSE contribution to secondary 
electron emission entirely. 

2 . Silver Bombarded by Ar , Ne, and Kr Ions 

The silver systems also indicate the possibility 
of the PSE mechanism being dependent on the crystal face 
being bombarded as well as the primary ion type. (See 
Figures 26 through 31.) 

In the Ar + - Ag system the PSE contribution to 
secondary electron emission is initially largest for the 
(100) face, followed by the (111) and (110) faces, respec- 
tively. In all three cases the PSE contribution decreases 
slightly at higher energies. 

In the Kr + - Ag system a different ordering of 
facial PSE contributions appear. Initially, the (111) 
face contributes the most to potential emission, followed 
by the (100) and (110) faces, in that order. At higher 
energies, the ordering changes, with (100) contributing 
the largest amount to PSE, followed by (111) and (110), 
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Figure 26. Ag bombarded by Ar + ions normally incident on the 
(100), (110) and (111) faces. 



73 




ENERGY (Kev) 

Figure 27. Ag bombarded by Ne + ions normally incident on the 
(100), (110) and (111) faces. 
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Figure 28. Ag bombarded by Kr + ions normally incident on the 
(100), (110) and (111) faces. 



PSE CONTRIBUTION 




ENERGY (KEV) 

Figure 29. PSE contribution to secondary electron emission 
for the Ar - Ag system. 
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Figure 30. PSE contribution to secondary electron emission 
for the Ne - Ag system. 
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Figure 31. PSE contribution to secondary electron emission 
. for the Kr - Ag system. 
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respectively. Once again, the PSE contribution tends to 
decrease at higher energies. 

In the Ne + - A g system an entirely different picture 
arises. For all energies, the (110) face contributes the 
most to PSE, followed by the (111) and (100) faces, 
respectively. However, even at ten Kev, the potential 
emission for the (110) face is still rising whereas, over 
the entire energy range, the PSE contribution from the 
(111) and (100) face remains fairly constant. These results 
would seem to indicate, once again, that not only is poten- 
tial emission dependent on primary ion type, but also on 
crystal face. 

In the Ar + - Ag system, the PSE contribution from 
the (111) face goes to a negative value at seven Kev. The 
reason for this behavior can be found by analyzing the 
Ar + - Ag SEE coefficient versus energy curves shown in 
Figure 26.- The SEE coefficient experimental data for the 
(111) face dips downward at seven Kev. This behavior is 
not demonstrated on the other faces and would tend to indi- 
cate that some experimental variation may be present in 
these experimental data. On the other hand, the interaction 
being modeled may be sufficiently complicated that the 
potential contribution cannot be separated simply by scaling 
the results. No other system investigated displays similar 
results . 
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3 . Molybdenum Bombarded by Ar , Ne , Kr and Xe Ions 



In Figures 32 through 43, Mo + initially seems to 
be the ion core model that best represents secondary electron 
emission from various target faces. This conclusion holds 
true for the Ne + - Mo and Kr + - Mo systems. However, for 
the Xe + - Mo system, the best ion core model is not Mo + but 
Mo ++ . This would seem to indicate that the primary ion 
type plays an important role in secondary electron emission. 

4 . Aluminum Bombarded by Ar , Kr , and Xe Ions 

In the case of aluminum bombarded by Ar + (Figures 44 
through 52) the Al + and Al ++ ion core models for the target 
atoms yielded results that were almost identical. This 
system is contrasted with the Kr + - A1 system which shows a 
distinct preference for Al ++ as the ion core model, with the 
Al + and Al +++ yielding similar and uniformly poor results. 
Similarly, the Xe + - A1 system shows the same preference 
for the Al ++ ion core model and demonstrates a poorer fit 
when the Al + and Al +++ models are used. When all cases are 
considered, it would seem that secondary electron emission 
is sensitive to the primary ion type, as has been demonstrated 
by previous results. 

C. ALKALI-HALIDE RESULTS 

Experimental results for the bombardment of alkali- 
halide single crystals by inert gas ions have been obtained 
from experiments by Baboux and Pedrix [21, 22]. Two cases 
are considered in this present investigation. These include 
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Figure 32. Effect of varying target atom ionization in the' 
HCM mo(jlel on secondary electron emission from 
the At - Mo system. 
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Figure 33. Effect of varying target atom ionization in the 
HCM mo^el on secondary electron emission from 
the Ar - Mo system. 
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Figure 34. Effect of varying target atom ionization in the 
HCy model on secondary electron emission in the 
, Ar - Mo system. 
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Figure 35. Effect on varying target atom ionization in the 
HCM mo^el on secondary electron emission from 
the Ne - Mo system. 
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Figure 36. Effect of varying target atom ionization in the 
HCM mocjel on secondary electron emission from 
* . the Ne - Mo system. 
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Figure 37. Effect of varying target atom ionization in the 
HCM mocjel on secondary electron emission from 
. the Ne - Mo system. 
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Figure 38. Effect of varying target atom ionization in 
the HCM mod^l on secondary electron emission 
from the Kr - Mo system. 
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Figure 39. Effect of varying target atom ionization in the 
HCM mo^Jel on secondary electron emission from 
. the Kr - Mo system. 
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Figure 40. Effect of varying target atom ionization in 
the HCM nod^l on secondary electron emission 
from tlu Kr - Mo system. 
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Figure 41. Effect of varying target atom ionization in 
the HCM mod^l on secondary electron emission 
from the Xe - Mo system. 
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Figure 42. Effect of varying target atom ionization in 
the HCM mod^l on secondary electron emission 
. from the Xe - Mo system. 
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Figure 43. Effect of varying target atom ionization in the 
HCM mo<Jel on secondary electron emission from 
. the Xe - Mo system. 



92 



T 



T 



T 



T 



T 




I I ! I I I I 1 I 

1 2 34 56789 10 
ENERGY(Kev) 

* » 

Figure 44. Effect of varying target atom ionization in the 
HCM mo^el on secondary electron emission from 
c # the Ar - A1 system. c 
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Figure 45. Effect of varying target atom ionization in the 

HCJjl model on secondary electron emission from the 
Ar - A1 system. 
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Figure 46. Effect of varying target atom ionization in the 
HCM mo^el on secondary electron emission from 
the Ar - Al system. 
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Figure 47. Effect of varying target atom ionization in the 
HCM mo^el on secondary electron emission from 
the Kr - A1 system. 
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Figure 48. Effect of varying target atom ionization in the 
HCM mocjel on secondary electron emission from 
the Kr - A1 system. 
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Figure 49. Effect of varying target atom ionization in the 
HCM mocjel on secondary electron emission from 
the Kr - A1 system. 
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Figure 50. Effect of varying target atom ionization in the' 
HCM mocjlel on secondary electron emission from 
the Xe - A1 system. 
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Figure 51. Effect of varying target atom ionization in the 
HCM mocjel on secondary electron emission from 
the Xe - A1 system. 
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Figure 52. Effect of varying target atom ionization in the 
HCM mocjel on secondary electron emission from 
the Xe - A1 system. 
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the bombardment of potassium chloride by argon and neon 
ions normally incident on the (100) , (110) and (111) faces 
of the single crystal. These two cases are discussed in 
detail below. 

1 . Potassium-Chloride Bombarded by Ar Ions 

The calculations of the SEE coefficient for alkali- 
halide single crystals present a theoretical problem 
that was not addressed when the metal target results 
were presented. Specifically, the alkali-halide crystals 
in question have lattice structures comprised of two dif- 
ferent atoms arranged in a similar fashion, but displaced 
in space by a distance (1/2, 1/2, 1/2). Because of this, 
the standard FCC distribution of impact parameters does 
not work in the alkali-halide case and must be replaced 
by a distribution representative of the normal impact 
of an ion with the target surface. The representative 
areas and distribution of impact parameters are chosen 
in a manner similar to that for the metal case, but are 
weighted to include preferential impact with the larger 
of the two ion cores in the lattice. The representative 
areas and distribution of impact parameters are shown in 
Figures 53 through 57. 

The results for KC1 bombarded by Ar + are shown in 
Figure 58. Experimentally, it was found that the SEE 
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□ Cl atoms 




K atoms 



Figure 54. Representative area KCI (110). 
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Figure 56. Potassium distribution of impact parameters 
in KC1 . 
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Figure 57. Chlorine distribution of impact parameters in 
* KC1 . 
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coefficient for the (100) face crossed that of the (110) 
face in the five to six Kev region. This behavior, as 
well as that for the (111) and (110) faces, is predicted 
by the HCM modified model. A PSE contribution is evident but 
has not been investigated. 

2. Potassium-Chloride Bombarded by Ne Ions 

In the case of KC1 bombarded by Ne + , the HCM 
modified model predicts a crossing of the (100) and (110) 
surfaces’ SEE coefficients as a function of energy at 
approximately five Kev. Such a behavior is not confirmed 
experimentally . 

In an attempt to explain this anomaly, two theories 
arose. First, such a behavior can be explained, if the KC1 
crystal tends to cleave on the (100) face with all the 
chlorine or potassium atoms up. Second, the crossing be- 
havior might be due to preferential sputtering by the primary 
ion of one of the two types of atoms comprising the (100) 
face, leaving a surface that contains either all potassiums 
or all chlorines. 

In the first case, the crystal structure of the (100) 
face should be the same for both Ar + and Ne + bombardment, 
as long as preferential sputtering does not occur. Hence, 

SEE from this face should be due to the bombardment of 
either pure chlorine or pure potassium, whichever is left 
in the "up” position. 

In the second case, however, preferential sputtering 
of chlorine atoms by Ar + is the case most likely to occur. 
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since both the ion and target atom are of comparable mass. 

In this situation, preferential sputtering would affect the 
Ar + - KC1 system more than it would the Ne + - KC1 system. 
Since the Ar + - KC1 system’s SEE coefficient as a function 
of energy is accurately predicted for all surfaces examined, 
it appears that sputtering is not directly responsible for 
the results that are obtained from either system. 

The effect of an all-potassium or all-chlorine 
approximation to the KC1 (100) surface on the SEE coefficient 
as a function of energy and primary ion type is shown graphi- 
cally in Figures 58 and 59. 
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Figure 59. SEE coef f icient + as a function of energy for KC1 
bombarded by Ne ions 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



This thesis has attempted to predict the behavior of the 
coefficient of ion electron emission for various primary ion- 
target combinations using the Harrison, Carlston and Magnuson 
theory of secondary electron emission. Three aspects of 
secondary electron emission were investigated, including 
angularly dependent secondary electron emission from the 
(100) face of a copper single crystal rotated about <110>, 
secondary electron emission from Ar , Ne, Kr and Xe ions nor- 
mally incident on the (100), (110) and (111) faces of Cu , Ag, 
Mo and A1 single crystals, and secondary electron emission 
from Ar and Ne ions normally incident on the (100) , (110) and 
(111) faces of KC1 . 

According to the HCM theory, angularly dependent results 
do not demonstrate the non-monotonic dependence on the angle 
of incidence found experimentally. Even so, the HCM theory 
does accurately model the dependence of the SEE coefficient 
of angle of incidence for inert gas ions normally incident on 
a polycrystalline metal surface. For this reason, it appears 
that the single collision model is not sophisticated enough 
and further modification must be made to it, if secondary 
electron emission from rotated, single-crystal faces is to be 
predicted. 

In the case of metal targets bombarded by inert gas 
ions, it was found that the SEE coefficient as a function of 
energy, predicted by the modified HCM movel, deviated from 
that experimentally obtained. This was contrasted with the 
results obtained using HCM’s original model which yielded 
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reasonable results in almost all cases investigated. Since 
the modifications presented here were refinements to the 
original HCM model, it was hoped the results obtained would 
more accurately represent the secondary electron mechanism. 

For this reason, the difference between the SEE coefficient 
obtained here and the experimental result was identified as 
the potential secondary electron emission component operating 
at higher energies. Such a hypothesis is substantiated by 
experiment [1], Analysis showed that the potential component 
in general rose to a peak and then decreased at higher energies. 
Furthermore, this component of secondary emission seemed to 
demonstrate behavior that could possibly be attributed to 
primary ion type as well as the face of the crystal being 
bombarded. The reason for this behavior was not positively 
identified, but the PSE process appears to be more complicated 
than was originally thought. 

In the case of potass ium- chloride bombarded by Ar ions, 
the HCM model fairly accurately predicted the SEE coefficient 
in the zero to 20 Kev energy range for all faces examined. 

On the other hand, potassium-chloride bombarded by Ne ions 
demonstrated a crossing of the SEE coefficient for the (100) 
and (110) faces as a function of energy which was not predicted 
experimentally. Two hypotheses were proposed that might pos- 
sibly explain this behavior. The first hypothesis proposes 
that KC1 cleaves with either all the potassium atoms or all 
the chlorine atoms up, essentially presenting a uniform 
layer of pure potassium or chlorine to the bombarding ion 
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beam. This hypothesis slightly improved the SEE 
coefficient for the Ne + KC1 (100) surface, but adversely 
affected the results obtained for the Ar + - KC1 system. 

The second hypothesis assumes that the (100) surface of KC1 
might demonstrate preferential sputtering. However, if 
this were the case, preferential sputtering should have been 
most evident in the Ar + - KC1 system, because of the compara 
ble mass of Ar and Cl. Since this system demonstrated the 
correct behavior without invoking the preferential sputter- 
ing assumption, it was surmised that preferential sputtering 
did not play an important role in the SEE process from the 
KC1 (100) surface. 

In all cases tested in the present investigation, the 
SEE coefficient as a function of energy tended to become a 
better approximation to the experimental results at higher 
energies. For this reason, it is recommended that primary 
ion-target atom combinations be run for primary ion energies 
between ten and 20 Kev. If in fact the potential component 
of secondary electron emission is real in this energy range, 
it should reach some constant value as primary ion energy 
is increased. Furthermore, because results obtained using 
the modified HCM model are dependent on the ion core configu 
ration of the target atoms, it is recommended that more 
work be done in testing the influence of ion-target atom 
core ionization on secondary electron emission. Typically, 
for secondary electron emission from elements such as Cu 
and Ag, which appear in column I of the periodic table, the 
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ion core electron configuration is uniquely determined by 
considerations of the atomic valence, but the effect of the 
electronic structure of the ion core upon KSE for Mo and 
A1 needs more investigation. 
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APPENDIX A 



COMPUTER PROGRAMS 



In this appendix, the parameters that appear in the 
computer programs are presented and explained. Variables 
that appear in the data input as constants are stated as 
such. A pictorial representation of how these data are 
entered into the programs is shown in Figures 12 - 16. 



A. HERMAN SKILLMAN RADIAL ELECTRON DENSITY PROGRAM 



CARD # 


NAME 


FORMAT 


DESCRIPTION 


1 


HEADER 


CARD - 


CONTAINS 80 CHARACTERS OF INFORMATION 
DESCRIBING THE CASE TO BE RUN 


2 


KEY 


(14) 


SEE Reference ( ) , 0 


2 


TOL 


(F8.6) 


" , 000001 


2 


THRESH 


(14) 


" , 521 


2 


MESH 


(14) 


" , 1 


2 


I PRATT 


(14) 


" , 40 


2 


MAX IT 


(14) 


" , 0 


2 


NOCOPY 


(14) 


" , 0 


2 


KUT 


(14) 


" , 0 


3-A* 


RU 2 ( I ) 


(IPE15 


.7, Normalized self-consistent poten- 



IP4E14.7) tial as a function of X. Taken 

directly from Reference 16. This 
is a "starter" for generating the 
wavefunctions for the atom under 
investigation. 

Atomic number of atom in question 
No. of shells of atom in question 
0 

Degree of ionization of atom in 
question 

n,l,m quantum no. for shells 
No. electrons in shell 
Trial eigenvalue for each shell, 
taken from Reference 16. 



B Z (F4.0) 

B NCORES (14) 

B NVALES (14) 

B ION (14) 



C NNLZ(I) (14) 

C WWNL(I) (F4.0) 

C EE (I ) (F8.4) 
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B. POTENTIAL PROGRAM 



CARD # 


NAME 


FORMAT 


DESCRIPTION 


1 


INUC 


(12) 


0 




1 


KSTEP 


(13) 


10 




1 


AZ 


(F3.0) 


Atomic no. of 


the first atom 






following the 


header card below 


1 


AI 


(F2.0) 


Ionization of 


the first atom 






following the 


header card. 


1 


BZ 


(F3. 0) 


Atomic no. of 


the second atom 






following the 


header card. 


1 


BI 


(F2.0) 


Ionization of 


the second atom 






following the 


header card. 


1 


DM 


(F2.0) 


5 




1 


RM 


(F2 . 0) 


5 




1 


RTEST 


(F2 . 0) 


5 




1 


DD 


(F5.1) 


.020 




1 


FRR 


(F5.1) 


1.20 




1 


DRS 


(F5.1) 


. 010 




1 


IPN 


(14) 


1 





2 HEADER CARD - CONTAINS 30 CHARACTERS OF INFORMATION 

DESCRIBING THE CASE TO BE RUN 



3 OUTPUT FROM ELECTRON DENSITY PROGRAM WITH FIRST 
CARD REMOVED CORRESPONDING TO THE ATOM DESCRIBED 
BY AZ AND AI ABOVE 

4 OUTPUT FROM ELECTRON DENSITY PROGRAM WITH FIRST 
CARD REMOVED CORRESPONDING TO THE ATOM DESCRIBED 
BY BZ AND BI ABOVE 



C. DISTRIBUTION OF IMPACT PARAMETERS PROGRAM 



CARD # 


NAME 


FORMAT 


DESCRIPTION 


1 


HEADER 


CARD - 80 CHARACTERS OF INFORMATION 
DESCRIBING THE 'CASE TO BE RUN 


2 


ICODE 


(13) 


1, if representative area is 
a rectangle 

2, if representative area is 
a triangle 


2 


A 


(F10.4) 


Height of representative area 
in Ro units. 


2 


B 


(F10 . 4) 


Length of representative area 
in Ro units. 


2 


INUMBER 


(13) 


No. of atoms in representative 
area . 
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C. DISTRIBUTION OF IMPACT PARAMETERS PROGRAM (Continued) 



CARD # NAME FORMAT DESCRIPTION 



3 


LARGE (I) 


(12) 


Logical array that indicates if 
any particular atom is larger 
than the other, i.e., Cl in 
KC1. 1, yes. 2, no. If yes 
the distribution is skewed to 
indicate preferential collision 
with the larger atom. 


4 


SCALEF 


(F10.4) 


1.0 


4 


RNA 


(F10.4) 


Ionic radius of smaller of two 
atoms . 


4 


RCL 


(F10.4) 


Ionic radius of larger of two 
atoms- -Both RNA and RCL can 
be set to anything if atoms 
are both the same size. 


5-7 


A5 ( I , J) 


(3F10.4) 


An array that contains the 
(x,y) vertices of a triangular 
representative area in Ro 
units. If the area is not tri- 
angular, they can be set to 
anything. 


8-11 


A6 (I , J) 


(3F10 . 4) 


An array that contains the (x,y 
position of atoms in the repre- 
sentative area in Ro units. I 
runs from 1 to INUMBR, J runs 
from 1 to 3. 


12 


PHI 


(F10.4) 


Angle in radians x f -axis makes 
with x-axis. 


12 


XR 


(F10.4) 


(x) position of (x’,y') origin. 


12 


YR 


(F10 . 4) 


(y) position of (x^y’) origin. 


12 


ALPHA 


(F10 . 4) 


Angle in radians the representa 
tive area is rotated through. 


12 


RMAX 


(F10 . 4) 


Maximum impact parameter 



expected. The no. 250 is 
usually sufficient. 
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D. NUMERICAL INTERPOLATION PROGRAM 



CARD # 


NAME 


FORMAT 


DESCRIPTION 


1 


NPOINT 


(14) 


No. points from the distribution 
of impact parameters program to 
be smoothed. 


1 


IFLAG 


(14) 


1, indicates written output only. 

2, indicates written output, 








card output and graphical 
output desired. 


2 


HEADER 


CARD - 46 


CHARACTERS TO BE OUTPUT WITH GRAPHS 


3 


HEADER 


CARD - 46 


CHARACTERS TO BE OUTPUT WITH GRAPHS 


4 


OUTPUT 


FROM DISTRIBUTION OF IMPACT PARAMETERS PRO- 




GRAM WITH ALL "NEW ATOM” CARDS REMOVED 


E. KSE 


PROGRAM 






CARD # 


NAME 


FORMAT 


DESCRIPTION 


1 


LCRYST 


(A4) 


Crystal type "FCC’V’BCC" 


1 


LTPE 


(12) 


1 


1 


BMAX ( I ) 


(3F6.4) 


Maximum impact parameter in Ro 
units for each face. 


1 


PUD (I ) 


(3F6.4) 


R^ for each face in Ro units. 


1 


ARO ( I ) 


(3F6. 4) 


No. atoms per representative 
area in Ro units, for each face. 


1 


LM(I) 


(313) 


Total number of points in each 
distribution. 


1 


KM 


(14) 


250 


1 


FACL 


(F7.4) 


(2) 1/2 /4 for FCC , (3) 1/2 /4 for 








BCC . 


2- A 


FLP(I) 


(2014) 


Distribution of impact parameters 
from Interpolation Program. 


B 


IH1 


80A1 


80 characters of information. 


C 


IH3 


19A4 


76 characters of information from 
first card of output from Potential 
Program. 


C 


KMR 


(14) 


250 


D 


APS I (I) 


10F8 . 5 


Wave functions from the Potential 
program, Atom A. 


D 


BPSI(I) 


10F8 . 5 


Wa' r e functions from the Potential 
program, Atom B. 
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E. KSE PROGRAM (Continued) 



CARD # 


NAME 


FORMAT 


E 


VCD 


5F12 . 5 


F 


EA(D 


8F8.4 


F+l 


I AH 


4A4 


F+l 


EA(D 


10F8.4 


F+2 


JFIT 


110 


G 


PSE(I) 


(F5.3) 



G 


FIT (J) 


(F5 . 4) 


H 


ZB 


F4.0 


H 


ZT 


F4 . 0 


H 


BMAS 


F8.4 


H 


TMAS 


F8.4 


II 


EOS 


F8 . 0 


II 


DEO 


F3. 0 


H 


KEOM 


F3 . 0 


H 


AO 


F7.3 


H 


LSKIP 


13 


H 


NM 


13 


H 


RES 


F3.0 


H 


DRE 


F3 . 0 


H 


REM 


F5 . 0 


H 


IH2 


4A4 



DESCRIPTION 



Potential between atoms A and 
B from the Potential Program. 

Spectroscopic ionization 
potentials, numbers 11 to 18. 
16 characters of accounting 
information . 

Spectroscopic ionization 
potentials, numbers 1 to 10. 
No. of values to which results 
are fit. 

YpSE experimental for curve 
being fit. 

y experimental for (111) face 
at 10 Kev. 

Atomic no. for atom A. 

Atomic No. for atom B. 

Atomic weight for atom A. 
Atomic weight for atom B. 
Energy in Kev at which results 
are to start. 

Energy step in Kev. 

Energy in Kev at which results 
are to end. 

Lattice constant in angstroms. 
1 

15 

50 

50 

5000 

Accounting data. 



^indicates that the actual card number varies from run to 
run and has been replaced by a letter to indicate the 
order in which the cards are arranged. 
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SLK2=SUK2+19.0*(VALUE+Y)+75.0*(SNLG(J+4)**2+SNLO( J+l)**2) 
+ 50. G*(SM_G{ J+2) **2+SNLO( J+3)**2 ) 

J = J + 5 
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THETA= i .57C7963 
R = StiRT ( A6 ( 1 ,2)**2+A6 l I , 1)**2) 

A6t It 1) =FX(R, THETA) 

A 6 ( I t 2) = FY(k,THETA) *C 0S ( ALPHA )-A6 ( I , 3 ) - S I N ( ALPHA) 
W F I Tt (6 ,t>5C) U A6 ( I, J ) , J=l,3) ,1 = 1 , INUMBR ) 



THIS SECTION CALCULATES THE ACTUAL DISTRIBUTION OF IMPACT PARAMETERS 
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SL6R0UT I N E CRTVAL EVALUATES THE ORTHOGONAL POLYNOMIALS AT 
SPECIFIED VALUE OF X 
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